#include <gmlib/number/mont256.h>
#include <gmlib/number/uint256.h>

#include <stdexcept>
#include <cstring>

using namespace number;

static std::uint8_t  P_data[32]      = {0xff,0xff,0xff,0xfe,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff};
static std::uint8_t  P_SUB2_data[32] = {0xff,0xff,0xff,0xfe,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xfd};
static std::uint8_t  R_data[32]      = {0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01};
static std::uint8_t  R_POW2_data[32] = {0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x02,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03};
static std::uint8_t  N_data[32]      = {0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01};

void test_mont256_dbl()
{
    uint256_t P, P_SUB2, R, R_POW2, N_;
    uint256_from_bytes(P, P_data);
    uint256_from_bytes(P_SUB2, P_SUB2_data);
    uint256_from_bytes(R, R_data);
    uint256_from_bytes(R_POW2, R_POW2_data);
    uint256_from_bytes(N_, N_data);

    Mont256CTX ctx;
    ctx.P = P, ctx.P_SUB2 = P_SUB2;
    ctx.R = R, ctx.R_POW2 = R_POW2;
    ctx.N_ = N_;

    mont256_t    multiplier, product;
    std::uint8_t product_data[32];

    std::uint8_t multiplier0[32] = {0xd2,0x3f,0x08,0x24,0x12,0x8b,0x2f,0x33,0x0c,0x5c,0x7f,0xd0,0xa6,0xa3,0xa4,0x50,0x65,0x13,0x27,0x0e,0x26,0x9e,0x0d,0x37,0xf2,0xa7,0x4d,0xe4,0x52,0xe6,0xb4,0x38};
    std::uint8_t product0[32]    = {0xa4,0x7e,0x10,0x49,0x25,0x16,0x5e,0x66,0x18,0xb8,0xff,0xa1,0x4d,0x47,0x48,0xa0,0xca,0x26,0x4e,0x1d,0x4d,0x3c,0x1a,0x6e,0xe5,0x4e,0x9b,0xc8,0xa5,0xcd,0x68,0x71};
    mont256_from_bytes(&ctx, multiplier, multiplier0);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product0, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier1[32] = {0x36,0xf6,0x75,0xcc,0x81,0xe7,0x4e,0xf5,0xe8,0xe2,0x5d,0x94,0x0e,0xd9,0x04,0x75,0x95,0x31,0x98,0x5d,0x5d,0x9d,0xc9,0xf8,0x18,0x18,0xe8,0x11,0x89,0x2f,0x90,0x2b};
    std::uint8_t product1[32]    = {0x6d,0xec,0xeb,0x99,0x03,0xce,0x9d,0xeb,0xd1,0xc4,0xbb,0x28,0x1d,0xb2,0x08,0xeb,0x2a,0x63,0x30,0xba,0xbb,0x3b,0x93,0xf0,0x30,0x31,0xd0,0x23,0x12,0x5f,0x20,0x56};
    mont256_from_bytes(&ctx, multiplier, multiplier1);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product1, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier2[32] = {0x8d,0x11,0x6e,0xce,0x17,0x38,0xf7,0xd9,0x3d,0x9c,0x17,0x24,0x11,0xe2,0x0b,0x8f,0x6b,0x0d,0x54,0x9b,0x6f,0x03,0x67,0x5a,0x16,0x00,0xa3,0x5a,0x09,0x99,0x50,0xd8};
    std::uint8_t product2[32]    = {0x1a,0x22,0xdd,0x9d,0x2e,0x71,0xef,0xb2,0x7b,0x38,0x2e,0x48,0x23,0xc4,0x17,0x1e,0xd6,0x1a,0xa9,0x37,0xde,0x06,0xce,0xb3,0x2c,0x01,0x46,0xb4,0x13,0x32,0xa1,0xb1};
    mont256_from_bytes(&ctx, multiplier, multiplier2);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product2, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier3[32] = {0xa1,0x70,0xb3,0x38,0x39,0x26,0x30,0x59,0xf2,0x8c,0x10,0x5d,0x1f,0xb1,0x7c,0x23,0x90,0xc1,0x92,0xcf,0xd3,0xac,0x94,0xaf,0x0f,0x21,0xdd,0xb6,0x6c,0xad,0x4a,0x26};
    std::uint8_t product3[32]    = {0x42,0xe1,0x66,0x71,0x72,0x4c,0x60,0xb3,0xe5,0x18,0x20,0xba,0x3f,0x62,0xf8,0x47,0x21,0x83,0x25,0xa0,0xa7,0x59,0x29,0x5d,0x1e,0x43,0xbb,0x6c,0xd9,0x5a,0x94,0x4d};
    mont256_from_bytes(&ctx, multiplier, multiplier3);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product3, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier4[32] = {0x0c,0xb1,0xe2,0x9c,0x65,0x8c,0xda,0x14,0x95,0xe6,0x0a,0xf5,0x93,0xbd,0x04,0xcf,0x0f,0xd6,0x30,0xf1,0xf2,0x9d,0x0d,0xa9,0x95,0x3f,0x48,0xf1,0xa0,0x9f,0x76,0xb5};
    std::uint8_t product4[32]    = {0x19,0x63,0xc5,0x38,0xcb,0x19,0xb4,0x29,0x2b,0xcc,0x15,0xeb,0x27,0x7a,0x09,0x9e,0x1f,0xac,0x61,0xe3,0xe5,0x3a,0x1b,0x53,0x2a,0x7e,0x91,0xe3,0x41,0x3e,0xed,0x6a};
    mont256_from_bytes(&ctx, multiplier, multiplier4);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product4, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier5[32] = {0x6b,0x4c,0xb2,0x42,0x4a,0x23,0xd5,0x96,0x22,0x17,0xbe,0xad,0xdb,0xc4,0x96,0xcb,0x8e,0x81,0x97,0x3e,0x0b,0xec,0xd7,0xb0,0x38,0x98,0xd1,0x90,0xf9,0xeb,0xda,0xcc};
    std::uint8_t product5[32]    = {0xd6,0x99,0x64,0x84,0x94,0x47,0xab,0x2c,0x44,0x2f,0x7d,0x5b,0xb7,0x89,0x2d,0x97,0x1d,0x03,0x2e,0x7c,0x17,0xd9,0xaf,0x60,0x71,0x31,0xa3,0x21,0xf3,0xd7,0xb5,0x98};
    mont256_from_bytes(&ctx, multiplier, multiplier5);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product5, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier6[32] = {0xae,0x97,0xba,0x94,0xd0,0xed,0xa8,0x2f,0x8f,0x6d,0x05,0x58,0x4e,0xf8,0xaa,0x38,0x92,0x27,0x66,0x58,0x1e,0x27,0xa1,0xc0,0x8a,0x6a,0x63,0xec,0x24,0xed,0xe6,0xa4};
    std::uint8_t product6[32]    = {0x5d,0x2f,0x75,0x2a,0xa1,0xdb,0x50,0x5f,0x1e,0xda,0x0a,0xb0,0x9d,0xf1,0x54,0x71,0x24,0x4e,0xcc,0xb1,0x3c,0x4f,0x43,0x80,0x14,0xd4,0xc7,0xd8,0x49,0xdb,0xcd,0x49};
    mont256_from_bytes(&ctx, multiplier, multiplier6);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product6, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier7[32] = {0x18,0xf1,0x35,0xd2,0x5f,0x55,0x72,0x03,0x30,0x18,0x50,0xc5,0xa3,0x8f,0xd5,0x47,0x92,0x3a,0x73,0x69,0x94,0xe3,0xbf,0x91,0x1a,0x61,0xdb,0xe2,0x2e,0x44,0x15,0x8b};
    std::uint8_t product7[32]    = {0x31,0xe2,0x6b,0xa4,0xbe,0xaa,0xe4,0x06,0x60,0x30,0xa1,0x8b,0x47,0x1f,0xaa,0x8f,0x24,0x74,0xe6,0xd3,0x29,0xc7,0x7f,0x22,0x34,0xc3,0xb7,0xc4,0x5c,0x88,0x2b,0x16};
    mont256_from_bytes(&ctx, multiplier, multiplier7);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product7, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier8[32] = {0x7f,0x15,0x05,0x24,0x34,0xb9,0xb5,0xdf,0x9e,0x77,0x69,0xb1,0x0f,0x42,0x05,0xb4,0x90,0x7a,0x70,0xc3,0x10,0x12,0xf0,0x37,0xb6,0x4c,0xe4,0x22,0x8c,0x38,0xfb,0x29};
    std::uint8_t product8[32]    = {0xfe,0x2a,0x0a,0x48,0x69,0x73,0x6b,0xbf,0x3c,0xee,0xd3,0x62,0x1e,0x84,0x0b,0x69,0x20,0xf4,0xe1,0x86,0x20,0x25,0xe0,0x6f,0x6c,0x99,0xc8,0x45,0x18,0x71,0xf6,0x52};
    mont256_from_bytes(&ctx, multiplier, multiplier8);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product8, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier9[32] = {0xec,0x66,0xa7,0x87,0x95,0xe7,0x61,0xd1,0x77,0x31,0xaf,0x10,0x50,0x6b,0xf2,0xef,0xc6,0xf8,0x77,0x18,0x6d,0x76,0xb0,0x7e,0x88,0x1e,0xd1,0x62,0xae,0x2e,0xb1,0x54};
    std::uint8_t product9[32]    = {0xd8,0xcd,0x4f,0x10,0x2b,0xce,0xc3,0xa2,0xee,0x63,0x5e,0x20,0xa0,0xd7,0xe5,0xdf,0x8d,0xf0,0xee,0x31,0xda,0xed,0x60,0xfc,0x10,0x3d,0xa2,0xc5,0x5c,0x5d,0x62,0xa9};
    mont256_from_bytes(&ctx, multiplier, multiplier9);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product9, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier10[32] = {0xc7,0xa2,0xea,0x20,0xb2,0xf1,0x4c,0x94,0x2e,0x05,0x31,0x9a,0xcb,0x5c,0x74,0x27,0x3f,0x98,0xe2,0x77,0x4c,0xbd,0x87,0xad,0x5c,0x90,0xa9,0x58,0x74,0x03,0xe4,0x30};
    std::uint8_t product10[32]    = {0x8f,0x45,0xd4,0x42,0x65,0xe2,0x99,0x28,0x5c,0x0a,0x63,0x35,0x96,0xb8,0xe8,0x4e,0x7f,0x31,0xc4,0xef,0x99,0x7b,0x0f,0x59,0xb9,0x21,0x52,0xb0,0xe8,0x07,0xc8,0x61};
    mont256_from_bytes(&ctx, multiplier, multiplier10);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product10, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier11[32] = {0x57,0xee,0x05,0xcd,0xe0,0x09,0x02,0xc7,0x7e,0xbf,0xf2,0x06,0x86,0x73,0x47,0x21,0x4c,0xdd,0x20,0x55,0x93,0x0d,0x6e,0xaf,0x14,0xf4,0x73,0x3f,0x3e,0x7d,0x1b,0xfb};
    std::uint8_t product11[32]    = {0xaf,0xdc,0x0b,0x9b,0xc0,0x12,0x05,0x8e,0xfd,0x7f,0xe4,0x0d,0x0c,0xe6,0x8e,0x42,0x99,0xba,0x40,0xab,0x26,0x1a,0xdd,0x5e,0x29,0xe8,0xe6,0x7e,0x7c,0xfa,0x37,0xf6};
    mont256_from_bytes(&ctx, multiplier, multiplier11);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product11, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier12[32] = {0x83,0x0e,0x07,0xbc,0x1e,0x39,0x8f,0x10,0x12,0xbd,0x4a,0xce,0xfa,0xec,0xbd,0x38,0x9b,0xe4,0xbc,0xfc,0x49,0xb6,0x4a,0x08,0x72,0xe6,0xcc,0x3a,0xba,0xbc,0xed,0x20};
    std::uint8_t product12[32]    = {0x06,0x1c,0x0f,0x79,0x3c,0x73,0x1e,0x20,0x25,0x7a,0x95,0x9d,0xf5,0xd9,0x7a,0x71,0x37,0xc9,0x79,0xf9,0x93,0x6c,0x94,0x0f,0xe5,0xcd,0x98,0x75,0x75,0x79,0xda,0x41};
    mont256_from_bytes(&ctx, multiplier, multiplier12);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product12, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier13[32] = {0x6b,0xf4,0x6c,0x69,0x7d,0x2c,0xaf,0x82,0xee,0xea,0xcb,0xe2,0x26,0xe8,0x75,0x55,0x57,0x90,0xf8,0x2e,0xc1,0xd3,0xfc,0xff,0x2a,0x3a,0xf4,0xd4,0x6b,0x0a,0x18,0xe8};
    std::uint8_t product13[32]    = {0xd7,0xe8,0xd8,0xd2,0xfa,0x59,0x5f,0x05,0xdd,0xd5,0x97,0xc4,0x4d,0xd0,0xea,0xaa,0xaf,0x21,0xf0,0x5d,0x83,0xa7,0xf9,0xfe,0x54,0x75,0xe9,0xa8,0xd6,0x14,0x31,0xd0};
    mont256_from_bytes(&ctx, multiplier, multiplier13);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product13, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier14[32] = {0xca,0x02,0x13,0x5e,0x92,0xb1,0xd3,0xf2,0x8e,0xde,0x0d,0x7a,0xc3,0xba,0xea,0x9e,0x13,0xde,0xef,0x86,0xab,0x10,0x31,0xd0,0xf6,0x46,0xe1,0xf4,0x0a,0x09,0x7c,0x97};
    std::uint8_t product14[32]    = {0x94,0x04,0x26,0xbe,0x25,0x63,0xa7,0xe5,0x1d,0xbc,0x1a,0xf5,0x87,0x75,0xd5,0x3c,0x27,0xbd,0xdf,0x0e,0x56,0x20,0x63,0xa0,0xec,0x8d,0xc3,0xe8,0x14,0x12,0xf9,0x2f};
    mont256_from_bytes(&ctx, multiplier, multiplier14);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product14, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier15[32] = {0x7f,0x26,0x14,0x4b,0x98,0x28,0x9f,0xcd,0x59,0xa5,0x4a,0x7b,0xb1,0xfe,0xe0,0x8f,0x57,0x12,0x42,0x42,0x50,0x51,0xc1,0xcc,0xd1,0x7f,0x9a,0xca,0xe0,0x1f,0x50,0x57};
    std::uint8_t product15[32]    = {0xfe,0x4c,0x28,0x97,0x30,0x51,0x3f,0x9a,0xb3,0x4a,0x94,0xf7,0x63,0xfd,0xc1,0x1e,0xae,0x24,0x84,0x84,0xa0,0xa3,0x83,0x99,0xa2,0xff,0x35,0x95,0xc0,0x3e,0xa0,0xae};
    mont256_from_bytes(&ctx, multiplier, multiplier15);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product15, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier16[32] = {0x45,0x1a,0xbd,0x81,0xf1,0xd6,0x9e,0xd6,0x17,0xf5,0xe8,0x37,0xd7,0x08,0x20,0xfe,0x11,0x9a,0x72,0xd1,0x74,0xc9,0xdf,0x6a,0xcc,0x01,0x1c,0xdd,0x94,0x74,0x03,0x1b};
    std::uint8_t product16[32]    = {0x8a,0x35,0x7b,0x03,0xe3,0xad,0x3d,0xac,0x2f,0xeb,0xd0,0x6f,0xae,0x10,0x41,0xfc,0x23,0x34,0xe5,0xa2,0xe9,0x93,0xbe,0xd5,0x98,0x02,0x39,0xbb,0x28,0xe8,0x06,0x36};
    mont256_from_bytes(&ctx, multiplier, multiplier16);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product16, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier17[32] = {0x4f,0x42,0x6d,0xcb,0xb3,0x94,0xfb,0x36,0xbb,0x2d,0x42,0x0f,0x0f,0x88,0x08,0x0b,0x10,0xa3,0xd6,0xb2,0xaa,0x05,0xe1,0x1a,0xb2,0x71,0x59,0x45,0x79,0x5e,0x82,0x29};
    std::uint8_t product17[32]    = {0x9e,0x84,0xdb,0x97,0x67,0x29,0xf6,0x6d,0x76,0x5a,0x84,0x1e,0x1f,0x10,0x10,0x16,0x21,0x47,0xad,0x65,0x54,0x0b,0xc2,0x35,0x64,0xe2,0xb2,0x8a,0xf2,0xbd,0x04,0x52};
    mont256_from_bytes(&ctx, multiplier, multiplier17);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product17, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier18[32] = {0xb7,0x74,0xeb,0x52,0x48,0xdb,0x40,0xaf,0x72,0x15,0x83,0x70,0xd2,0x69,0xa9,0xa5,0xae,0x65,0x8f,0x33,0xfe,0x3b,0x89,0x0b,0x93,0xf4,0x48,0xb3,0xa5,0xaa,0x3c,0x81};
    std::uint8_t product18[32]    = {0x6e,0xe9,0xd6,0xa5,0x91,0xb6,0x81,0x5e,0xe4,0x2b,0x06,0xe1,0xa4,0xd3,0x53,0x4b,0x5c,0xcb,0x1e,0x68,0xfc,0x77,0x12,0x16,0x27,0xe8,0x91,0x67,0x4b,0x54,0x79,0x03};
    mont256_from_bytes(&ctx, multiplier, multiplier18);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product18, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier19[32] = {0x5a,0xff,0xb2,0x29,0x76,0x31,0xa9,0x92,0xf0,0xce,0x58,0x35,0x05,0xc6,0xaf,0x07,0x58,0xd5,0x56,0x3d,0xab,0x2c,0xd3,0x1e,0xe3,0x15,0x12,0x88,0x62,0xc3,0x3a,0x4f};
    std::uint8_t product19[32]    = {0xb5,0xff,0x64,0x52,0xec,0x63,0x53,0x25,0xe1,0x9c,0xb0,0x6a,0x0b,0x8d,0x5e,0x0e,0xb1,0xaa,0xac,0x7b,0x56,0x59,0xa6,0x3d,0xc6,0x2a,0x25,0x10,0xc5,0x86,0x74,0x9e};
    mont256_from_bytes(&ctx, multiplier, multiplier19);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product19, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier20[32] = {0x49,0x95,0x23,0x99,0xc4,0xaa,0xea,0xc1,0x37,0xdc,0x76,0xfb,0x0f,0x17,0xa3,0x00,0x7e,0x62,0xaa,0x0a,0x1d,0xf9,0xfd,0x78,0x9c,0x65,0x39,0x38,0x2b,0x05,0x37,0xe6};
    std::uint8_t product20[32]    = {0x93,0x2a,0x47,0x33,0x89,0x55,0xd5,0x82,0x6f,0xb8,0xed,0xf6,0x1e,0x2f,0x46,0x00,0xfc,0xc5,0x54,0x14,0x3b,0xf3,0xfa,0xf1,0x38,0xca,0x72,0x70,0x56,0x0a,0x6f,0xcc};
    mont256_from_bytes(&ctx, multiplier, multiplier20);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product20, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier21[32] = {0x7f,0x1b,0x10,0x3c,0xdf,0x15,0x82,0xb0,0xea,0xb4,0x77,0xd2,0x64,0x15,0x47,0x9c,0x65,0xdc,0x9f,0x50,0x3f,0x63,0xaf,0x83,0xbd,0x05,0x61,0xe6,0x21,0x1c,0x70,0xcf};
    std::uint8_t product21[32]    = {0xfe,0x36,0x20,0x79,0xbe,0x2b,0x05,0x61,0xd5,0x68,0xef,0xa4,0xc8,0x2a,0x8f,0x38,0xcb,0xb9,0x3e,0xa0,0x7e,0xc7,0x5f,0x07,0x7a,0x0a,0xc3,0xcc,0x42,0x38,0xe1,0x9e};
    mont256_from_bytes(&ctx, multiplier, multiplier21);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product21, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier22[32] = {0x23,0x0d,0x97,0x7e,0xe2,0x25,0x71,0x59,0x47,0x20,0x77,0x1f,0x8c,0xa8,0x18,0x11,0x66,0xd2,0x28,0x76,0x72,0xfd,0xf2,0x02,0x2a,0x96,0xfb,0x1a,0x14,0xa0,0xf9,0xe7};
    std::uint8_t product22[32]    = {0x46,0x1b,0x2e,0xfd,0xc4,0x4a,0xe2,0xb2,0x8e,0x40,0xee,0x3f,0x19,0x50,0x30,0x22,0xcd,0xa4,0x50,0xec,0xe5,0xfb,0xe4,0x04,0x55,0x2d,0xf6,0x34,0x29,0x41,0xf3,0xce};
    mont256_from_bytes(&ctx, multiplier, multiplier22);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product22, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier23[32] = {0xfc,0x89,0x1b,0x4a,0x6a,0x50,0xdf,0x4d,0xb4,0xd6,0x6a,0x3a,0x47,0x46,0x9a,0x4d,0x8c,0xdb,0x30,0x5f,0xdd,0x2e,0x16,0x09,0x6e,0x36,0xaa,0xb0,0xd1,0xbc,0x52,0xd9};
    std::uint8_t product23[32]    = {0xf9,0x12,0x36,0x95,0xd4,0xa1,0xbe,0x9b,0x69,0xac,0xd4,0x74,0x8e,0x8d,0x34,0x9b,0x19,0xb6,0x60,0xc0,0xba,0x5c,0x2c,0x11,0xdc,0x6d,0x55,0x61,0xa3,0x78,0xa5,0xb3};
    mont256_from_bytes(&ctx, multiplier, multiplier23);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product23, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier24[32] = {0x15,0x3e,0x7c,0x2a,0x26,0xa2,0xc0,0xbd,0x3b,0x12,0x87,0xff,0xf5,0x2d,0xdf,0x5d,0x61,0x64,0x99,0xc9,0xe2,0x5a,0x76,0x05,0xae,0xc6,0xf0,0x24,0x5b,0xd8,0x6d,0x40};
    std::uint8_t product24[32]    = {0x2a,0x7c,0xf8,0x54,0x4d,0x45,0x81,0x7a,0x76,0x25,0x0f,0xff,0xea,0x5b,0xbe,0xba,0xc2,0xc9,0x33,0x93,0xc4,0xb4,0xec,0x0b,0x5d,0x8d,0xe0,0x48,0xb7,0xb0,0xda,0x80};
    mont256_from_bytes(&ctx, multiplier, multiplier24);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product24, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier25[32] = {0xd4,0xc2,0x8c,0x2e,0x7c,0x26,0x84,0x7f,0x03,0x16,0x90,0x9e,0x3b,0xbb,0xe9,0xea,0xa8,0x94,0x8c,0x89,0x3b,0x61,0x86,0x76,0x26,0xbb,0x7d,0xbd,0x2d,0x1c,0x9a,0xf0};
    std::uint8_t product25[32]    = {0xa9,0x85,0x18,0x5d,0xf8,0x4d,0x08,0xfe,0x06,0x2d,0x21,0x3c,0x77,0x77,0xd3,0xd5,0x51,0x29,0x19,0x13,0x76,0xc3,0x0c,0xeb,0x4d,0x76,0xfb,0x7a,0x5a,0x39,0x35,0xe1};
    mont256_from_bytes(&ctx, multiplier, multiplier25);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product25, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier26[32] = {0x88,0xda,0xf4,0x01,0x6b,0x40,0x13,0xef,0x25,0x4b,0x0c,0x4e,0x01,0x0c,0x47,0x59,0x48,0x2c,0x9c,0xbc,0x43,0x43,0x5c,0xc5,0x2e,0xae,0x05,0xcf,0x96,0xd0,0xcc,0x5f};
    std::uint8_t product26[32]    = {0x11,0xb5,0xe8,0x03,0xd6,0x80,0x27,0xde,0x4a,0x96,0x18,0x9c,0x02,0x18,0x8e,0xb2,0x90,0x59,0x39,0x79,0x86,0x86,0xb9,0x89,0x5d,0x5c,0x0b,0x9f,0x2d,0xa1,0x98,0xbf};
    mont256_from_bytes(&ctx, multiplier, multiplier26);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product26, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier27[32] = {0xdb,0xf4,0xa8,0xb2,0xb0,0xc4,0x31,0x2d,0x20,0x20,0x36,0x26,0xf3,0xfe,0x39,0xc0,0x51,0x90,0x88,0xf5,0x90,0xfb,0xbd,0x11,0x9c,0x1c,0xaa,0xf7,0x5e,0x87,0x66,0xed};
    std::uint8_t product27[32]    = {0xb7,0xe9,0x51,0x66,0x61,0x88,0x62,0x5a,0x40,0x40,0x6c,0x4d,0xe7,0xfc,0x73,0x80,0xa3,0x21,0x11,0xec,0x21,0xf7,0x7a,0x22,0x38,0x39,0x55,0xee,0xbd,0x0e,0xcd,0xdb};
    mont256_from_bytes(&ctx, multiplier, multiplier27);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product27, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier28[32] = {0x74,0xe6,0x9a,0x5d,0x0d,0xd2,0x7a,0x65,0xbd,0x62,0x88,0x81,0xad,0x1b,0x72,0xdb,0xa7,0xab,0xe1,0xc2,0x9e,0x1a,0x8e,0xf4,0xf3,0x41,0xe0,0x7a,0x83,0xf7,0x3f,0x16};
    std::uint8_t product28[32]    = {0xe9,0xcd,0x34,0xba,0x1b,0xa4,0xf4,0xcb,0x7a,0xc5,0x11,0x03,0x5a,0x36,0xe5,0xb7,0x4f,0x57,0xc3,0x85,0x3c,0x35,0x1d,0xe9,0xe6,0x83,0xc0,0xf5,0x07,0xee,0x7e,0x2c};
    mont256_from_bytes(&ctx, multiplier, multiplier28);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product28, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier29[32] = {0x8f,0x2c,0x6e,0xc8,0xcc,0x41,0x69,0xa3,0xae,0x3a,0x2b,0x7f,0xdf,0xe0,0x18,0x93,0xf3,0xae,0xd0,0xb6,0xc7,0xac,0x14,0x91,0xde,0xf8,0x83,0x34,0xe6,0x47,0xcb,0x8f};
    std::uint8_t product29[32]    = {0x1e,0x58,0xdd,0x92,0x98,0x82,0xd3,0x47,0x5c,0x74,0x56,0xff,0xbf,0xc0,0x31,0x27,0xe7,0x5d,0xa1,0x6e,0x8f,0x58,0x29,0x22,0xbd,0xf1,0x06,0x69,0xcc,0x8f,0x97,0x1f};
    mont256_from_bytes(&ctx, multiplier, multiplier29);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product29, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier30[32] = {0x66,0x83,0x68,0x86,0xa2,0x60,0xcd,0x0b,0x7b,0x45,0x14,0x5c,0x1a,0x81,0x68,0x2c,0x64,0xe5,0x0c,0xad,0x66,0x23,0x7a,0x04,0x65,0xe7,0xe4,0x23,0x64,0x72,0xf1,0xa3};
    std::uint8_t product30[32]    = {0xcd,0x06,0xd1,0x0d,0x44,0xc1,0x9a,0x16,0xf6,0x8a,0x28,0xb8,0x35,0x02,0xd0,0x58,0xc9,0xca,0x19,0x5a,0xcc,0x46,0xf4,0x08,0xcb,0xcf,0xc8,0x46,0xc8,0xe5,0xe3,0x46};
    mont256_from_bytes(&ctx, multiplier, multiplier30);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product30, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier31[32] = {0x1c,0x24,0x42,0xf9,0x29,0x8c,0xb3,0xa5,0x70,0xcc,0xec,0x31,0x35,0x71,0x81,0x0a,0xfc,0x13,0x2d,0x0d,0x11,0x3d,0xb1,0x7d,0x30,0xcb,0xc9,0x7d,0x0f,0xef,0x79,0x28};
    std::uint8_t product31[32]    = {0x38,0x48,0x85,0xf2,0x53,0x19,0x67,0x4a,0xe1,0x99,0xd8,0x62,0x6a,0xe3,0x02,0x15,0xf8,0x26,0x5a,0x1a,0x22,0x7b,0x62,0xfa,0x61,0x97,0x92,0xfa,0x1f,0xde,0xf2,0x50};
    mont256_from_bytes(&ctx, multiplier, multiplier31);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product31, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier32[32] = {0x89,0x5f,0xd7,0xb3,0x26,0xb9,0x4c,0x7f,0x91,0x18,0xbb,0x16,0x00,0x0f,0x49,0xc8,0x1a,0x35,0x8c,0xa0,0x0d,0x75,0x98,0x5d,0x99,0xc9,0x43,0x09,0x57,0x0d,0xc1,0x95};
    std::uint8_t product32[32]    = {0x12,0xbf,0xaf,0x67,0x4d,0x72,0x98,0xff,0x22,0x31,0x76,0x2c,0x00,0x1e,0x93,0x90,0x34,0x6b,0x19,0x41,0x1a,0xeb,0x30,0xba,0x33,0x92,0x86,0x12,0xae,0x1b,0x83,0x2b};
    mont256_from_bytes(&ctx, multiplier, multiplier32);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product32, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier33[32] = {0x35,0x3c,0x63,0x1c,0xdf,0xd4,0x3f,0x37,0x12,0x00,0x33,0x9d,0x06,0x87,0x39,0xfa,0x9d,0x1d,0xe2,0xa0,0x5d,0x15,0x8a,0x2f,0xf2,0xee,0x4e,0x45,0x19,0xf9,0x91,0x9c};
    std::uint8_t product33[32]    = {0x6a,0x78,0xc6,0x39,0xbf,0xa8,0x7e,0x6e,0x24,0x00,0x67,0x3a,0x0d,0x0e,0x73,0xf5,0x3a,0x3b,0xc5,0x40,0xba,0x2b,0x14,0x5f,0xe5,0xdc,0x9c,0x8a,0x33,0xf3,0x23,0x38};
    mont256_from_bytes(&ctx, multiplier, multiplier33);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product33, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier34[32] = {0x9a,0x2e,0xf8,0x0f,0x58,0xee,0x85,0x71,0xf4,0x99,0x8d,0x7c,0x40,0x93,0xf6,0xde,0xa2,0x68,0xaa,0x87,0x26,0x07,0x67,0x9d,0x60,0x50,0x91,0x4a,0x9d,0x33,0xa0,0x1c};
    std::uint8_t product34[32]    = {0x34,0x5d,0xf0,0x1f,0xb1,0xdd,0x0a,0xe3,0xe9,0x33,0x1a,0xf8,0x81,0x27,0xed,0xbd,0x44,0xd1,0x55,0x0f,0x4c,0x0e,0xcf,0x39,0xc0,0xa1,0x22,0x95,0x3a,0x67,0x40,0x39};
    mont256_from_bytes(&ctx, multiplier, multiplier34);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product34, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier35[32] = {0xfa,0x52,0x9b,0xa3,0xfe,0x3b,0xfa,0xda,0x7c,0xf2,0x07,0x24,0xd9,0x53,0xee,0x26,0x1d,0x87,0xce,0xc3,0x1f,0x72,0x96,0xab,0x79,0x61,0xfd,0x92,0x5d,0x39,0xd0,0xa8};
    std::uint8_t product35[32]    = {0xf4,0xa5,0x37,0x48,0xfc,0x77,0xf5,0xb4,0xf9,0xe4,0x0e,0x49,0xb2,0xa7,0xdc,0x4c,0x3b,0x0f,0x9d,0x87,0x3e,0xe5,0x2d,0x55,0xf2,0xc3,0xfb,0x24,0xba,0x73,0xa1,0x51};
    mont256_from_bytes(&ctx, multiplier, multiplier35);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product35, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier36[32] = {0xbf,0xea,0xa1,0x55,0x1a,0x28,0xf7,0xb3,0x24,0xe4,0xe2,0x5a,0x15,0xfc,0x89,0x9e,0x4f,0xd5,0x8d,0xbe,0x7b,0xdc,0x96,0x8b,0x7a,0xfb,0x2c,0x68,0x77,0x4b,0x15,0xd7};
    std::uint8_t product36[32]    = {0x7f,0xd5,0x42,0xab,0x34,0x51,0xef,0x66,0x49,0xc9,0xc4,0xb4,0x2b,0xf9,0x13,0x3c,0x9f,0xab,0x1b,0x7d,0xf7,0xb9,0x2d,0x15,0xf5,0xf6,0x58,0xd0,0xee,0x96,0x2b,0xaf};
    mont256_from_bytes(&ctx, multiplier, multiplier36);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product36, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier37[32] = {0x84,0x2e,0x7f,0xc2,0x29,0x54,0x0a,0x6e,0xb1,0x2a,0xa1,0xf6,0xd4,0x2f,0xdd,0xbb,0x7a,0x86,0xf7,0xa2,0x43,0xc7,0x1b,0x9a,0xbd,0x87,0xa8,0x65,0x57,0xb6,0xfb,0x7e};
    std::uint8_t product37[32]    = {0x08,0x5c,0xff,0x85,0x52,0xa8,0x14,0xdd,0x62,0x55,0x43,0xed,0xa8,0x5f,0xbb,0x76,0xf5,0x0d,0xef,0x45,0x87,0x8e,0x37,0x34,0x7b,0x0f,0x50,0xca,0xaf,0x6d,0xf6,0xfd};
    mont256_from_bytes(&ctx, multiplier, multiplier37);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product37, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier38[32] = {0xb0,0xa8,0x44,0xe5,0x25,0x87,0xbe,0x6b,0x5c,0x9b,0xcf,0x35,0x87,0x3b,0xe0,0x78,0xf3,0xb7,0xa5,0x0d,0xf3,0x73,0xca,0x53,0x34,0x88,0xf8,0x76,0x05,0xe9,0x99,0xf3};
    std::uint8_t product38[32]    = {0x61,0x50,0x89,0xcb,0x4b,0x0f,0x7c,0xd6,0xb9,0x37,0x9e,0x6b,0x0e,0x77,0xc0,0xf1,0xe7,0x6f,0x4a,0x1c,0xe6,0xe7,0x94,0xa5,0x69,0x11,0xf0,0xec,0x0b,0xd3,0x33,0xe7};
    mont256_from_bytes(&ctx, multiplier, multiplier38);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product38, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier39[32] = {0xa4,0x96,0x36,0xa2,0xfa,0x7f,0x0e,0xab,0x4c,0x4f,0x9b,0x06,0x87,0x32,0x2e,0x25,0xc2,0x15,0xa8,0x2a,0x06,0xec,0x41,0xad,0xea,0x05,0x75,0x43,0x8b,0x0d,0x59,0x0b};
    std::uint8_t product39[32]    = {0x49,0x2c,0x6d,0x46,0xf4,0xfe,0x1d,0x56,0x98,0x9f,0x36,0x0d,0x0e,0x64,0x5c,0x4b,0x84,0x2b,0x50,0x55,0x0d,0xd8,0x83,0x5a,0xd4,0x0a,0xea,0x87,0x16,0x1a,0xb2,0x17};
    mont256_from_bytes(&ctx, multiplier, multiplier39);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product39, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier40[32] = {0xe8,0x83,0xa1,0xd4,0x5d,0xe0,0x09,0x97,0x84,0xb5,0xa8,0x18,0x42,0xd8,0x72,0x08,0xd8,0x6f,0x40,0xf6,0xb2,0x39,0xf3,0xc7,0x17,0x4c,0x77,0xa2,0xdd,0x02,0xde,0x92};
    std::uint8_t product40[32]    = {0xd1,0x07,0x43,0xa9,0xbb,0xc0,0x13,0x2f,0x09,0x6b,0x50,0x30,0x85,0xb0,0xe4,0x11,0xb0,0xde,0x81,0xee,0x64,0x73,0xe7,0x8d,0x2e,0x98,0xef,0x45,0xba,0x05,0xbd,0x25};
    mont256_from_bytes(&ctx, multiplier, multiplier40);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product40, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier41[32] = {0x80,0xb0,0xc0,0x8b,0xc7,0x70,0x24,0x20,0x8a,0xa4,0x24,0x8c,0x88,0x57,0xf9,0xa4,0x39,0x08,0xf2,0x27,0xc5,0x9d,0xb9,0x16,0x5b,0x0e,0xe7,0x6f,0x2a,0xc3,0x44,0x46};
    std::uint8_t product41[32]    = {0x01,0x61,0x81,0x18,0x8e,0xe0,0x48,0x41,0x15,0x48,0x49,0x19,0x10,0xaf,0xf3,0x48,0x72,0x11,0xe4,0x50,0x8b,0x3b,0x72,0x2b,0xb6,0x1d,0xce,0xde,0x55,0x86,0x88,0x8d};
    mont256_from_bytes(&ctx, multiplier, multiplier41);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product41, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier42[32] = {0xc2,0x21,0x6b,0x02,0xfc,0x24,0x1d,0x0b,0xc9,0xd4,0x88,0xb1,0xcf,0xbf,0x33,0x60,0x9c,0xfc,0x86,0x52,0x39,0x19,0x42,0x42,0xa2,0xed,0xdb,0xbd,0x54,0x64,0xec,0xc2};
    std::uint8_t product42[32]    = {0x84,0x42,0xd6,0x06,0xf8,0x48,0x3a,0x17,0x93,0xa9,0x11,0x63,0x9f,0x7e,0x66,0xc1,0x39,0xf9,0x0c,0xa5,0x72,0x32,0x84,0x84,0x45,0xdb,0xb7,0x7a,0xa8,0xc9,0xd9,0x85};
    mont256_from_bytes(&ctx, multiplier, multiplier42);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product42, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier43[32] = {0xcd,0xa6,0xc6,0xfd,0xbd,0x68,0x51,0x67,0x66,0x93,0x40,0x36,0xd1,0x7e,0x44,0x97,0x3d,0x48,0x82,0xa5,0xce,0x5b,0x2a,0x92,0x31,0xf5,0x17,0x07,0xda,0x45,0xe1,0x8a};
    std::uint8_t product43[32]    = {0x9b,0x4d,0x8d,0xfc,0x7a,0xd0,0xa2,0xce,0xcd,0x26,0x80,0x6d,0xa2,0xfc,0x89,0x2e,0x7a,0x91,0x05,0x4c,0x9c,0xb6,0x55,0x23,0x63,0xea,0x2e,0x0f,0xb4,0x8b,0xc3,0x15};
    mont256_from_bytes(&ctx, multiplier, multiplier43);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product43, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier44[32] = {0xfd,0x56,0xa9,0x26,0x07,0x6b,0x3e,0x36,0xbb,0x23,0x13,0xf5,0x5b,0x06,0x25,0x8e,0x7e,0x26,0xf3,0x6a,0x84,0x83,0xf8,0xb8,0x33,0x2d,0xd3,0x31,0x3a,0x0b,0x99,0x65};
    std::uint8_t product44[32]    = {0xfa,0xad,0x52,0x4d,0x0e,0xd6,0x7c,0x6d,0x76,0x46,0x27,0xea,0xb6,0x0c,0x4b,0x1c,0xfc,0x4d,0xe6,0xd6,0x09,0x07,0xf1,0x6f,0x66,0x5b,0xa6,0x62,0x74,0x17,0x32,0xcb};
    mont256_from_bytes(&ctx, multiplier, multiplier44);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product44, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier45[32] = {0x9a,0xea,0x64,0x29,0xb1,0x49,0x1e,0x24,0x31,0x92,0xb7,0x04,0x42,0x59,0x40,0x52,0x78,0xe4,0xb9,0x8d,0x47,0x87,0xf9,0x3b,0xca,0x44,0xeb,0x86,0x07,0x26,0xe2,0x5c};
    std::uint8_t product45[32]    = {0x35,0xd4,0xc8,0x54,0x62,0x92,0x3c,0x48,0x63,0x25,0x6e,0x08,0x84,0xb2,0x80,0xa4,0xf1,0xc9,0x73,0x1b,0x8f,0x0f,0xf2,0x76,0x94,0x89,0xd7,0x0c,0x0e,0x4d,0xc4,0xb9};
    mont256_from_bytes(&ctx, multiplier, multiplier45);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product45, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier46[32] = {0x59,0x7a,0x1e,0xcf,0xfc,0xf0,0x0f,0xec,0xb9,0x1e,0xe9,0xe5,0xef,0xe0,0x9f,0x07,0xce,0xfe,0x2a,0x1f,0x72,0x7d,0x83,0x49,0x58,0x22,0xcb,0x77,0xf4,0xde,0x2c,0x08};
    std::uint8_t product46[32]    = {0xb2,0xf4,0x3d,0x9f,0xf9,0xe0,0x1f,0xd9,0x72,0x3d,0xd3,0xcb,0xdf,0xc1,0x3e,0x0f,0x9d,0xfc,0x54,0x3e,0xe4,0xfb,0x06,0x92,0xb0,0x45,0x96,0xef,0xe9,0xbc,0x58,0x10};
    mont256_from_bytes(&ctx, multiplier, multiplier46);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product46, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier47[32] = {0x78,0x57,0x29,0x76,0x3a,0x12,0x91,0x7c,0x1a,0x26,0xf8,0x89,0x38,0x70,0x38,0x00,0x14,0x9e,0x25,0x9b,0x5d,0x58,0xc7,0x05,0xf9,0x79,0xd0,0x4a,0xf4,0x7a,0xeb,0xdd};
    std::uint8_t product47[32]    = {0xf0,0xae,0x52,0xec,0x74,0x25,0x22,0xf8,0x34,0x4d,0xf1,0x12,0x70,0xe0,0x70,0x00,0x29,0x3c,0x4b,0x36,0xba,0xb1,0x8e,0x0b,0xf2,0xf3,0xa0,0x95,0xe8,0xf5,0xd7,0xba};
    mont256_from_bytes(&ctx, multiplier, multiplier47);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product47, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier48[32] = {0x9c,0x3a,0x23,0xcd,0xe6,0x7a,0x9b,0x75,0xfc,0x39,0x47,0x24,0x9f,0xc2,0xd0,0xa1,0x7b,0x8f,0x2a,0xb5,0x34,0x51,0xd0,0x13,0x56,0x75,0xf6,0xad,0x32,0x5b,0x55,0xdd};
    std::uint8_t product48[32]    = {0x38,0x74,0x47,0x9c,0xcc,0xf5,0x36,0xeb,0xf8,0x72,0x8e,0x49,0x3f,0x85,0xa1,0x42,0xf7,0x1e,0x55,0x6b,0x68,0xa3,0xa0,0x25,0xac,0xeb,0xed,0x5a,0x64,0xb6,0xab,0xbb};
    mont256_from_bytes(&ctx, multiplier, multiplier48);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product48, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier49[32] = {0xa4,0xa4,0x5e,0xff,0xcc,0xb5,0x73,0xd9,0x58,0x10,0xd6,0x0e,0xa7,0x29,0x91,0xb9,0xe8,0xc1,0x47,0x43,0x7a,0xbe,0xc5,0x39,0x00,0x7d,0x10,0x34,0xd7,0x26,0xc8,0x6b};
    std::uint8_t product49[32]    = {0x49,0x48,0xbe,0x00,0x99,0x6a,0xe7,0xb2,0xb0,0x21,0xac,0x1d,0x4e,0x53,0x23,0x73,0xd1,0x82,0x8e,0x87,0xf5,0x7d,0x8a,0x71,0x00,0xfa,0x20,0x69,0xae,0x4d,0x90,0xd7};
    mont256_from_bytes(&ctx, multiplier, multiplier49);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product49, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier50[32] = {0xb6,0x24,0x67,0x71,0xc8,0x45,0x00,0x70,0x63,0x77,0x14,0x07,0xe8,0xe7,0x27,0x89,0x1e,0xb2,0x01,0x09,0xa9,0x1c,0x24,0x39,0xd5,0xab,0x8b,0x4d,0x15,0xb4,0x0a,0xeb};
    std::uint8_t product50[32]    = {0x6c,0x48,0xce,0xe4,0x90,0x8a,0x00,0xe0,0xc6,0xee,0x28,0x0f,0xd1,0xce,0x4f,0x12,0x3d,0x64,0x02,0x14,0x52,0x38,0x48,0x72,0xab,0x57,0x16,0x9a,0x2b,0x68,0x15,0xd7};
    mont256_from_bytes(&ctx, multiplier, multiplier50);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product50, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier51[32] = {0xa2,0xc6,0x8e,0x45,0xca,0x04,0xc7,0x9f,0x6f,0x15,0xb6,0xad,0x2d,0xb3,0x99,0x7f,0xe3,0x96,0x39,0xbe,0x7a,0x60,0x5a,0x91,0x33,0x06,0x98,0xa1,0xc0,0x09,0x34,0x92};
    std::uint8_t product51[32]    = {0x45,0x8d,0x1c,0x8c,0x94,0x09,0x8f,0x3e,0xde,0x2b,0x6d,0x5a,0x5b,0x67,0x32,0xff,0xc7,0x2c,0x73,0x7d,0xf4,0xc0,0xb5,0x21,0x66,0x0d,0x31,0x43,0x80,0x12,0x69,0x25};
    mont256_from_bytes(&ctx, multiplier, multiplier51);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product51, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier52[32] = {0x76,0x91,0xb0,0x6f,0x65,0x55,0xab,0xfe,0xb8,0xc9,0x81,0x7a,0xf8,0xbe,0x88,0x31,0xf2,0x37,0xe4,0x5a,0xcd,0x02,0xc5,0xe1,0x16,0x35,0x3d,0x03,0x55,0x1f,0xd8,0xf9};
    std::uint8_t product52[32]    = {0xed,0x23,0x60,0xde,0xca,0xab,0x57,0xfd,0x71,0x93,0x02,0xf5,0xf1,0x7d,0x10,0x63,0xe4,0x6f,0xc8,0xb5,0x9a,0x05,0x8b,0xc2,0x2c,0x6a,0x7a,0x06,0xaa,0x3f,0xb1,0xf2};
    mont256_from_bytes(&ctx, multiplier, multiplier52);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product52, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier53[32] = {0xfe,0x3c,0x9c,0x8f,0x2b,0x85,0x5c,0x1f,0x28,0xaa,0xca,0x51,0xb9,0x8c,0x67,0xc2,0x15,0xbd,0x44,0x8f,0xf2,0x61,0x49,0xed,0xbe,0x4c,0x5c,0xe6,0x66,0xc1,0x49,0x4e};
    std::uint8_t product53[32]    = {0xfc,0x79,0x39,0x1f,0x57,0x0a,0xb8,0x3e,0x51,0x55,0x94,0xa3,0x73,0x18,0xcf,0x84,0x2b,0x7a,0x89,0x20,0xe4,0xc2,0x93,0xda,0x7c,0x98,0xb9,0xcc,0xcd,0x82,0x92,0x9d};
    mont256_from_bytes(&ctx, multiplier, multiplier53);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product53, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier54[32] = {0xa7,0xe6,0x52,0x9b,0xce,0x76,0xe9,0xf4,0x77,0x21,0x6e,0x9e,0xe7,0xa4,0x63,0x09,0x97,0x3f,0x79,0x86,0x26,0xb1,0xcf,0xfc,0x07,0x0d,0x71,0x09,0x20,0x85,0x96,0x34};
    std::uint8_t product54[32]    = {0x4f,0xcc,0xa5,0x38,0x9c,0xed,0xd3,0xe8,0xee,0x42,0xdd,0x3d,0xcf,0x48,0xc6,0x13,0x2e,0x7e,0xf3,0x0d,0x4d,0x63,0x9f,0xf7,0x0e,0x1a,0xe2,0x12,0x41,0x0b,0x2c,0x69};
    mont256_from_bytes(&ctx, multiplier, multiplier54);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product54, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier55[32] = {0xef,0xfd,0xde,0xea,0xa8,0x42,0xbc,0x19,0x79,0x6f,0x74,0xad,0xfa,0xf5,0x54,0x96,0x98,0x8a,0xf3,0xfb,0xd3,0x96,0x30,0xd6,0x9c,0x90,0x11,0xef,0x25,0x6b,0xad,0xf9};
    std::uint8_t product55[32]    = {0xdf,0xfb,0xbd,0xd6,0x50,0x85,0x78,0x32,0xf2,0xde,0xe9,0x5b,0xf5,0xea,0xa9,0x2d,0x31,0x15,0xe7,0xf8,0xa7,0x2c,0x61,0xac,0x39,0x20,0x23,0xde,0x4a,0xd7,0x5b,0xf3};
    mont256_from_bytes(&ctx, multiplier, multiplier55);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product55, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier56[32] = {0xcc,0xa2,0xa9,0x2b,0x03,0xa5,0x6c,0xc1,0x05,0x7a,0x40,0xb2,0x21,0x88,0x28,0x7e,0x8c,0x5c,0x71,0x5f,0x8c,0x74,0xfc,0x1e,0x27,0xe9,0xe0,0x6f,0x59,0xb4,0x4e,0x92};
    std::uint8_t product56[32]    = {0x99,0x45,0x52,0x57,0x07,0x4a,0xd9,0x82,0x0a,0xf4,0x81,0x64,0x43,0x10,0x50,0xfd,0x18,0xb8,0xe2,0xc0,0x18,0xe9,0xf8,0x3b,0x4f,0xd3,0xc0,0xde,0xb3,0x68,0x9d,0x25};
    mont256_from_bytes(&ctx, multiplier, multiplier56);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product56, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier57[32] = {0x23,0xa5,0xef,0x88,0xef,0x02,0x09,0x0b,0xbf,0xde,0xfc,0x15,0x86,0xce,0x03,0xf9,0x1a,0x4f,0x44,0xf9,0xa6,0x51,0x14,0x45,0xb9,0xf3,0x63,0x5c,0xf8,0x8c,0x42,0x2b};
    std::uint8_t product57[32]    = {0x47,0x4b,0xdf,0x11,0xde,0x04,0x12,0x17,0x7f,0xbd,0xf8,0x2b,0x0d,0x9c,0x07,0xf2,0x34,0x9e,0x89,0xf3,0x4c,0xa2,0x28,0x8b,0x73,0xe6,0xc6,0xb9,0xf1,0x18,0x84,0x56};
    mont256_from_bytes(&ctx, multiplier, multiplier57);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product57, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier58[32] = {0x07,0x2a,0x98,0xd2,0x36,0x06,0xde,0xfc,0xdf,0xb8,0x5c,0x0d,0xd3,0x7e,0xe9,0x15,0x31,0xde,0xc4,0xf4,0xdf,0x2a,0x8b,0x79,0xfc,0x8e,0x80,0xb3,0x6f,0x0e,0x22,0x89};
    std::uint8_t product58[32]    = {0x0e,0x55,0x31,0xa4,0x6c,0x0d,0xbd,0xf9,0xbf,0x70,0xb8,0x1b,0xa6,0xfd,0xd2,0x2a,0x63,0xbd,0x89,0xe9,0xbe,0x55,0x16,0xf3,0xf9,0x1d,0x01,0x66,0xde,0x1c,0x45,0x12};
    mont256_from_bytes(&ctx, multiplier, multiplier58);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product58, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier59[32] = {0x53,0x74,0x09,0x02,0x96,0x20,0xbf,0x0d,0xc3,0x80,0x84,0xa0,0x3d,0x93,0xfd,0x4c,0x80,0x4c,0x25,0xd6,0x4a,0xff,0xdc,0xd1,0x36,0x78,0xbc,0x8d,0x40,0x78,0x3f,0x0a};
    std::uint8_t product59[32]    = {0xa6,0xe8,0x12,0x05,0x2c,0x41,0x7e,0x1b,0x87,0x01,0x09,0x40,0x7b,0x27,0xfa,0x99,0x00,0x98,0x4b,0xac,0x95,0xff,0xb9,0xa2,0x6c,0xf1,0x79,0x1a,0x80,0xf0,0x7e,0x14};
    mont256_from_bytes(&ctx, multiplier, multiplier59);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product59, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier60[32] = {0xbd,0x6b,0x88,0x1a,0xe8,0xf6,0xe0,0xbd,0x0f,0x97,0x70,0x44,0x21,0x8e,0x0b,0x7b,0xd5,0x8d,0xcd,0xb4,0x6b,0x44,0x68,0x06,0x8b,0x5a,0xb3,0xee,0x42,0x65,0xbb,0x31};
    std::uint8_t product60[32]    = {0x7a,0xd7,0x10,0x36,0xd1,0xed,0xc1,0x7a,0x1f,0x2e,0xe0,0x88,0x43,0x1c,0x16,0xf7,0xab,0x1b,0x9b,0x69,0xd6,0x88,0xd0,0x0c,0x16,0xb5,0x67,0xdc,0x84,0xcb,0x76,0x63};
    mont256_from_bytes(&ctx, multiplier, multiplier60);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product60, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier61[32] = {0x84,0x4a,0x70,0x34,0xe7,0x7f,0xfe,0x48,0xd0,0xa6,0xec,0x17,0x95,0x56,0x58,0x5e,0xa9,0x97,0xf3,0x51,0x75,0x4a,0x09,0xcd,0xe5,0xcf,0xed,0xfa,0x5a,0x91,0x96,0xf0};
    std::uint8_t product61[32]    = {0x08,0x94,0xe0,0x6a,0xce,0xff,0xfc,0x91,0xa1,0x4d,0xd8,0x2f,0x2a,0xac,0xb0,0xbd,0x53,0x2f,0xe6,0xa3,0xea,0x94,0x13,0x9a,0xcb,0x9f,0xdb,0xf4,0xb5,0x23,0x2d,0xe1};
    mont256_from_bytes(&ctx, multiplier, multiplier61);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product61, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier62[32] = {0x26,0xde,0xbf,0xdb,0x88,0x25,0xae,0x56,0x21,0x79,0xb3,0x7d,0x80,0x6c,0x10,0xb5,0xe0,0xcf,0xab,0x4c,0xea,0xef,0xc4,0xd2,0xd3,0xbf,0x6d,0x01,0x6b,0xae,0x4b,0x5b};
    std::uint8_t product62[32]    = {0x4d,0xbd,0x7f,0xb7,0x10,0x4b,0x5c,0xac,0x42,0xf3,0x66,0xfb,0x00,0xd8,0x21,0x6b,0xc1,0x9f,0x56,0x99,0xd5,0xdf,0x89,0xa5,0xa7,0x7e,0xda,0x02,0xd7,0x5c,0x96,0xb6};
    mont256_from_bytes(&ctx, multiplier, multiplier62);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product62, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier63[32] = {0x9b,0xca,0x3c,0xb7,0x2e,0xe0,0x28,0x9d,0xc6,0xc9,0x1b,0x92,0x70,0xac,0x06,0xac,0xdf,0x70,0x30,0x17,0x04,0xc9,0xd7,0x8d,0x82,0xb3,0x35,0x99,0x86,0x04,0x87,0x19};
    std::uint8_t product63[32]    = {0x37,0x94,0x79,0x6f,0x5d,0xc0,0x51,0x3b,0x8d,0x92,0x37,0x24,0xe1,0x58,0x0d,0x59,0xbe,0xe0,0x60,0x2f,0x09,0x93,0xaf,0x1a,0x05,0x66,0x6b,0x33,0x0c,0x09,0x0e,0x33};
    mont256_from_bytes(&ctx, multiplier, multiplier63);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product63, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier64[32] = {0x9e,0x7d,0x6b,0x37,0x79,0x36,0xd5,0x36,0x24,0x3d,0x35,0x70,0x2c,0x1e,0xea,0x1f,0x26,0x59,0x74,0xa7,0xcc,0x96,0x6f,0x46,0xc6,0xaa,0x7d,0x55,0x01,0x01,0xb8,0x11};
    std::uint8_t product64[32]    = {0x3c,0xfa,0xd6,0x6f,0xf2,0x6d,0xaa,0x6c,0x48,0x7a,0x6a,0xe0,0x58,0x3d,0xd4,0x3e,0x4c,0xb2,0xe9,0x50,0x99,0x2c,0xde,0x8c,0x8d,0x54,0xfa,0xaa,0x02,0x03,0x70,0x23};
    mont256_from_bytes(&ctx, multiplier, multiplier64);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product64, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier65[32] = {0x87,0xdd,0xae,0xb7,0x84,0xb2,0x80,0x54,0xae,0xad,0x44,0xb0,0x53,0x73,0x90,0xe5,0x0f,0xcf,0x31,0xca,0x8e,0x75,0x2f,0xdf,0x1e,0xce,0x61,0x5d,0xb9,0xa6,0x44,0x2e};
    std::uint8_t product65[32]    = {0x0f,0xbb,0x5d,0x70,0x09,0x65,0x00,0xa9,0x5d,0x5a,0x89,0x60,0xa6,0xe7,0x21,0xca,0x1f,0x9e,0x63,0x96,0x1c,0xea,0x5f,0xbd,0x3d,0x9c,0xc2,0xbb,0x73,0x4c,0x88,0x5d};
    mont256_from_bytes(&ctx, multiplier, multiplier65);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product65, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier66[32] = {0x0e,0x8b,0xec,0x94,0x8f,0x6f,0x91,0x5f,0xe2,0x1b,0x37,0xca,0x1b,0x29,0xfc,0x99,0xc6,0xc8,0x0e,0x2b,0xc8,0xc6,0x14,0xb2,0x7b,0x84,0x44,0xd1,0x8e,0x31,0x70,0x41};
    std::uint8_t product66[32]    = {0x1d,0x17,0xd9,0x29,0x1e,0xdf,0x22,0xbf,0xc4,0x36,0x6f,0x94,0x36,0x53,0xf9,0x33,0x8d,0x90,0x1c,0x57,0x91,0x8c,0x29,0x64,0xf7,0x08,0x89,0xa3,0x1c,0x62,0xe0,0x82};
    mont256_from_bytes(&ctx, multiplier, multiplier66);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product66, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier67[32] = {0x73,0xc1,0xcd,0x2c,0x81,0xf9,0x8b,0x52,0x19,0x05,0xd5,0x91,0xc5,0xb2,0xe7,0x5a,0x0a,0xcd,0x8b,0xe1,0x46,0xe4,0x09,0x90,0x30,0xf9,0x70,0x58,0x3f,0x9d,0x52,0xf9};
    std::uint8_t product67[32]    = {0xe7,0x83,0x9a,0x59,0x03,0xf3,0x16,0xa4,0x32,0x0b,0xab,0x23,0x8b,0x65,0xce,0xb4,0x15,0x9b,0x17,0xc2,0x8d,0xc8,0x13,0x20,0x61,0xf2,0xe0,0xb0,0x7f,0x3a,0xa5,0xf2};
    mont256_from_bytes(&ctx, multiplier, multiplier67);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product67, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier68[32] = {0x53,0x5b,0x6a,0x43,0x71,0x78,0xba,0x0a,0x10,0x38,0xf0,0xb5,0xe9,0x98,0xd0,0xee,0xe4,0xdd,0xf9,0xb9,0xc2,0x8e,0xe9,0x07,0x07,0x22,0x35,0xc2,0x8f,0xcd,0x7f,0x40};
    std::uint8_t product68[32]    = {0xa6,0xb6,0xd4,0x86,0xe2,0xf1,0x74,0x14,0x20,0x71,0xe1,0x6b,0xd3,0x31,0xa1,0xdd,0xc9,0xbb,0xf3,0x73,0x85,0x1d,0xd2,0x0e,0x0e,0x44,0x6b,0x85,0x1f,0x9a,0xfe,0x80};
    mont256_from_bytes(&ctx, multiplier, multiplier68);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product68, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier69[32] = {0x46,0xf5,0xa1,0xb4,0xb1,0x56,0xd1,0xad,0x33,0x0c,0x16,0xa3,0x83,0x1d,0x03,0xbf,0x9b,0x2b,0xd6,0xc0,0x81,0x6b,0xee,0x06,0xf9,0x2e,0x23,0x39,0x9c,0xce,0xa0,0x98};
    std::uint8_t product69[32]    = {0x8d,0xeb,0x43,0x69,0x62,0xad,0xa3,0x5a,0x66,0x18,0x2d,0x47,0x06,0x3a,0x07,0x7f,0x36,0x57,0xad,0x81,0x02,0xd7,0xdc,0x0d,0xf2,0x5c,0x46,0x73,0x39,0x9d,0x41,0x30};
    mont256_from_bytes(&ctx, multiplier, multiplier69);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product69, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier70[32] = {0x3f,0x66,0x5e,0xde,0xf1,0x06,0x37,0xce,0x81,0xfc,0x06,0x9e,0x7a,0x60,0x96,0x83,0xce,0xaf,0x49,0x15,0x88,0x85,0x64,0xe8,0x82,0x16,0x85,0x8f,0x73,0xcc,0xef,0x03};
    std::uint8_t product70[32]    = {0x7e,0xcc,0xbd,0xbd,0xe2,0x0c,0x6f,0x9d,0x03,0xf8,0x0d,0x3c,0xf4,0xc1,0x2d,0x07,0x9d,0x5e,0x92,0x2b,0x11,0x0a,0xc9,0xd1,0x04,0x2d,0x0b,0x1e,0xe7,0x99,0xde,0x06};
    mont256_from_bytes(&ctx, multiplier, multiplier70);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product70, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier71[32] = {0xec,0x3b,0x96,0x05,0x42,0x74,0xa3,0xeb,0xed,0x84,0xe9,0x1e,0xf1,0x32,0xbf,0x2d,0xe0,0x40,0x01,0x5c,0xe0,0x64,0xa1,0x14,0x85,0xf1,0x11,0x5b,0xb2,0xff,0xf1,0x7b};
    std::uint8_t product71[32]    = {0xd8,0x77,0x2c,0x0b,0x84,0xe9,0x47,0xd7,0xdb,0x09,0xd2,0x3d,0xe2,0x65,0x7e,0x5b,0xc0,0x80,0x02,0xba,0xc0,0xc9,0x42,0x28,0x0b,0xe2,0x22,0xb7,0x65,0xff,0xe2,0xf7};
    mont256_from_bytes(&ctx, multiplier, multiplier71);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product71, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier72[32] = {0x6a,0xa8,0xb9,0xe0,0x23,0x1b,0x3e,0x14,0x72,0x91,0x35,0xbd,0xd7,0x0a,0x39,0xd1,0x33,0xdc,0xd7,0x7f,0xf1,0x79,0xf2,0xd2,0xe4,0x8b,0x96,0x62,0x8f,0x3c,0x4b,0xe3};
    std::uint8_t product72[32]    = {0xd5,0x51,0x73,0xc0,0x46,0x36,0x7c,0x28,0xe5,0x22,0x6b,0x7b,0xae,0x14,0x73,0xa2,0x67,0xb9,0xae,0xff,0xe2,0xf3,0xe5,0xa5,0xc9,0x17,0x2c,0xc5,0x1e,0x78,0x97,0xc6};
    mont256_from_bytes(&ctx, multiplier, multiplier72);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product72, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier73[32] = {0x6d,0xa7,0x9a,0x87,0x3d,0x9a,0x80,0x79,0xab,0xd0,0xd7,0xfb,0x12,0x92,0x61,0x85,0x50,0xe4,0x0d,0x54,0x71,0x2e,0xa6,0xb3,0x64,0x71,0xfd,0xe4,0x1f,0x22,0x9d,0xd0};
    std::uint8_t product73[32]    = {0xdb,0x4f,0x35,0x0e,0x7b,0x35,0x00,0xf3,0x57,0xa1,0xaf,0xf6,0x25,0x24,0xc3,0x0a,0xa1,0xc8,0x1a,0xa8,0xe2,0x5d,0x4d,0x66,0xc8,0xe3,0xfb,0xc8,0x3e,0x45,0x3b,0xa0};
    mont256_from_bytes(&ctx, multiplier, multiplier73);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product73, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier74[32] = {0xc6,0xe5,0x0d,0xf2,0xe5,0xa3,0x86,0x3e,0x1f,0x52,0x52,0x65,0xc8,0xb0,0x07,0xee,0x4d,0x82,0xfe,0xac,0xab,0x62,0x86,0xcd,0x36,0x72,0xd6,0xae,0x12,0xb8,0x0a,0xed};
    std::uint8_t product74[32]    = {0x8d,0xca,0x1b,0xe6,0xcb,0x47,0x0c,0x7c,0x3e,0xa4,0xa4,0xcb,0x91,0x60,0x0f,0xdc,0x9b,0x05,0xfd,0x5a,0x56,0xc5,0x0d,0x99,0x6c,0xe5,0xad,0x5c,0x25,0x70,0x15,0xdb};
    mont256_from_bytes(&ctx, multiplier, multiplier74);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product74, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier75[32] = {0x40,0xcb,0xac,0xd0,0x24,0x9a,0x45,0x84,0x5d,0xbe,0x30,0x23,0xa9,0x06,0x92,0x2f,0xa4,0xb9,0xa9,0xc4,0xb7,0x53,0xa1,0xee,0xf0,0x83,0x60,0x85,0x27,0x89,0xd0,0x59};
    std::uint8_t product75[32]    = {0x81,0x97,0x59,0xa0,0x49,0x34,0x8b,0x08,0xbb,0x7c,0x60,0x47,0x52,0x0d,0x24,0x5f,0x49,0x73,0x53,0x89,0x6e,0xa7,0x43,0xdd,0xe1,0x06,0xc1,0x0a,0x4f,0x13,0xa0,0xb2};
    mont256_from_bytes(&ctx, multiplier, multiplier75);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product75, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier76[32] = {0x18,0x18,0x9a,0xf4,0xf3,0xd7,0x4f,0x82,0xbf,0x26,0x8e,0xa0,0x38,0x36,0xe8,0x65,0x77,0xbd,0x89,0x1f,0xf7,0xb1,0x03,0xdf,0x23,0x23,0x1e,0x1e,0xe2,0x01,0x55,0x22};
    std::uint8_t product76[32]    = {0x30,0x31,0x35,0xe9,0xe7,0xae,0x9f,0x05,0x7e,0x4d,0x1d,0x40,0x70,0x6d,0xd0,0xca,0xef,0x7b,0x12,0x3f,0xef,0x62,0x07,0xbe,0x46,0x46,0x3c,0x3d,0xc4,0x02,0xaa,0x44};
    mont256_from_bytes(&ctx, multiplier, multiplier76);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product76, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier77[32] = {0x39,0x45,0x33,0x6b,0xd5,0x1b,0x18,0x15,0xaa,0xf7,0x19,0xf3,0xfd,0x68,0x37,0x3b,0x29,0xac,0xf1,0xa5,0x7c,0xbd,0x1f,0x5a,0xe2,0x8a,0xf6,0x04,0x65,0xf4,0x29,0x86};
    std::uint8_t product77[32]    = {0x72,0x8a,0x66,0xd7,0xaa,0x36,0x30,0x2b,0x55,0xee,0x33,0xe7,0xfa,0xd0,0x6e,0x76,0x53,0x59,0xe3,0x4a,0xf9,0x7a,0x3e,0xb5,0xc5,0x15,0xec,0x08,0xcb,0xe8,0x53,0x0c};
    mont256_from_bytes(&ctx, multiplier, multiplier77);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product77, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier78[32] = {0x6b,0xd8,0xc6,0x76,0x56,0xd0,0x50,0xcd,0x67,0x60,0x13,0x67,0x83,0xfe,0xb1,0x7b,0xfe,0x7b,0x8a,0xe4,0x6e,0x78,0x36,0xa4,0xb4,0xd1,0x9e,0xc1,0x29,0x55,0xd6,0xf0};
    std::uint8_t product78[32]    = {0xd7,0xb1,0x8c,0xec,0xad,0xa0,0xa1,0x9a,0xce,0xc0,0x26,0xcf,0x07,0xfd,0x62,0xf7,0xfc,0xf7,0x15,0xc8,0xdc,0xf0,0x6d,0x49,0x69,0xa3,0x3d,0x82,0x52,0xab,0xad,0xe0};
    mont256_from_bytes(&ctx, multiplier, multiplier78);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product78, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier79[32] = {0x56,0x85,0xd6,0x24,0x04,0xfc,0xd5,0x55,0x5d,0xaf,0x10,0x6d,0xb8,0xde,0xe0,0x81,0x17,0x9a,0x07,0x1e,0x51,0x8a,0xe4,0x52,0x5b,0x4b,0x1b,0x75,0x32,0x1c,0x52,0x96};
    std::uint8_t product79[32]    = {0xad,0x0b,0xac,0x48,0x09,0xf9,0xaa,0xaa,0xbb,0x5e,0x20,0xdb,0x71,0xbd,0xc1,0x02,0x2f,0x34,0x0e,0x3c,0xa3,0x15,0xc8,0xa4,0xb6,0x96,0x36,0xea,0x64,0x38,0xa5,0x2c};
    mont256_from_bytes(&ctx, multiplier, multiplier79);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product79, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier80[32] = {0x84,0x76,0x8b,0x8c,0x54,0xdd,0x0b,0xa5,0x62,0x64,0x67,0xba,0x04,0xa1,0x05,0x47,0xb4,0x01,0xba,0x85,0x70,0xc1,0xdc,0xa1,0x75,0x6b,0x72,0x89,0x8d,0xd6,0x3c,0xb9};
    std::uint8_t product80[32]    = {0x08,0xed,0x17,0x19,0xa9,0xba,0x17,0x4a,0xc4,0xc8,0xcf,0x74,0x09,0x42,0x0a,0x8f,0x68,0x03,0x75,0x0b,0xe1,0x83,0xb9,0x41,0xea,0xd6,0xe5,0x13,0x1b,0xac,0x79,0x73};
    mont256_from_bytes(&ctx, multiplier, multiplier80);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product80, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier81[32] = {0xeb,0x25,0xf8,0xa1,0xfc,0x2e,0x6a,0x59,0x1c,0xe3,0xbc,0x0c,0x10,0x75,0x5c,0x97,0xf5,0xf5,0x54,0xed,0x83,0x23,0x9e,0xf5,0x4b,0xa2,0xe1,0x61,0x9f,0xb9,0xaf,0x50};
    std::uint8_t product81[32]    = {0xd6,0x4b,0xf1,0x44,0xf8,0x5c,0xd4,0xb2,0x39,0xc7,0x78,0x18,0x20,0xea,0xb9,0x2f,0xeb,0xea,0xa9,0xdc,0x06,0x47,0x3d,0xe9,0x97,0x45,0xc2,0xc3,0x3f,0x73,0x5e,0xa1};
    mont256_from_bytes(&ctx, multiplier, multiplier81);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product81, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier82[32] = {0x45,0x9c,0x94,0x5c,0x43,0xfc,0x05,0x27,0x15,0x85,0x0a,0x03,0x1a,0xd2,0xd5,0xf1,0xe0,0x5b,0x3e,0x13,0xf8,0xc1,0x10,0xfb,0x3a,0x82,0x81,0x59,0xc9,0xd2,0x29,0x50};
    std::uint8_t product82[32]    = {0x8b,0x39,0x28,0xb8,0x87,0xf8,0x0a,0x4e,0x2b,0x0a,0x14,0x06,0x35,0xa5,0xab,0xe3,0xc0,0xb6,0x7c,0x27,0xf1,0x82,0x21,0xf6,0x75,0x05,0x02,0xb3,0x93,0xa4,0x52,0xa0};
    mont256_from_bytes(&ctx, multiplier, multiplier82);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product82, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier83[32] = {0xd1,0xdc,0xec,0x53,0x21,0x2a,0x8d,0x9b,0xc1,0x7a,0x92,0x62,0x45,0x3b,0xf4,0x91,0x2e,0x7a,0x26,0xe9,0xc7,0x6c,0x60,0x3f,0xe7,0xe8,0xf9,0xf6,0x0a,0x22,0x73,0x85};
    std::uint8_t product83[32]    = {0xa3,0xb9,0xd8,0xa7,0x42,0x55,0x1b,0x37,0x82,0xf5,0x24,0xc4,0x8a,0x77,0xe9,0x22,0x5c,0xf4,0x4d,0xd4,0x8e,0xd8,0xc0,0x7e,0xcf,0xd1,0xf3,0xec,0x14,0x44,0xe7,0x0b};
    mont256_from_bytes(&ctx, multiplier, multiplier83);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product83, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier84[32] = {0x67,0xec,0x32,0x6a,0x42,0x34,0x33,0x54,0xf2,0x2d,0x28,0x82,0xd1,0xa8,0x9b,0x37,0xad,0x0c,0x9b,0xb6,0xe9,0x52,0x6a,0x69,0xd9,0x7e,0x96,0x7b,0x6c,0x18,0xd9,0x82};
    std::uint8_t product84[32]    = {0xcf,0xd8,0x64,0xd4,0x84,0x68,0x66,0xa9,0xe4,0x5a,0x51,0x05,0xa3,0x51,0x36,0x6f,0x5a,0x19,0x37,0x6d,0xd2,0xa4,0xd4,0xd3,0xb2,0xfd,0x2c,0xf6,0xd8,0x31,0xb3,0x04};
    mont256_from_bytes(&ctx, multiplier, multiplier84);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product84, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier85[32] = {0x53,0xb9,0x73,0x77,0xb3,0x4e,0x8e,0xce,0x7e,0x9e,0xe5,0x1d,0x92,0x12,0x82,0x4c,0x83,0xc8,0xcb,0x28,0xeb,0x4e,0xd2,0xe3,0x89,0x5e,0x8b,0x6b,0x26,0x3c,0xfa,0x5e};
    std::uint8_t product85[32]    = {0xa7,0x72,0xe6,0xef,0x66,0x9d,0x1d,0x9c,0xfd,0x3d,0xca,0x3b,0x24,0x25,0x04,0x99,0x07,0x91,0x96,0x51,0xd6,0x9d,0xa5,0xc7,0x12,0xbd,0x16,0xd6,0x4c,0x79,0xf4,0xbc};
    mont256_from_bytes(&ctx, multiplier, multiplier85);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product85, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier86[32] = {0xe5,0x31,0x69,0x60,0x6c,0xe1,0x93,0xc2,0x2e,0xef,0xa2,0x79,0xb0,0x2e,0x3d,0x8d,0xcc,0xb1,0xc5,0x1d,0x0e,0xba,0x0e,0xa8,0x47,0x70,0xa0,0x87,0x16,0xe6,0xfe,0xc3};
    std::uint8_t product86[32]    = {0xca,0x62,0xd2,0xc1,0xd9,0xc3,0x27,0x84,0x5d,0xdf,0x44,0xf3,0x60,0x5c,0x7b,0x1b,0x99,0x63,0x8a,0x3b,0x1d,0x74,0x1d,0x4f,0x8e,0xe1,0x41,0x0e,0x2d,0xcd,0xfd,0x87};
    mont256_from_bytes(&ctx, multiplier, multiplier86);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product86, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier87[32] = {0x42,0xb3,0x87,0x55,0xcd,0x37,0x88,0x0e,0x16,0xac,0x41,0x91,0xa2,0x6a,0xa0,0xae,0x04,0x4f,0x15,0x74,0xf0,0x37,0xaf,0xc6,0x44,0xd8,0x2a,0x53,0x12,0x89,0xba,0xfa};
    std::uint8_t product87[32]    = {0x85,0x67,0x0e,0xab,0x9a,0x6f,0x10,0x1c,0x2d,0x58,0x83,0x23,0x44,0xd5,0x41,0x5c,0x08,0x9e,0x2a,0xe9,0xe0,0x6f,0x5f,0x8c,0x89,0xb0,0x54,0xa6,0x25,0x13,0x75,0xf4};
    mont256_from_bytes(&ctx, multiplier, multiplier87);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product87, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier88[32] = {0x1f,0x26,0x42,0xaa,0xdc,0xde,0xd2,0x04,0x43,0xb3,0x0f,0x66,0x11,0x0e,0x2c,0xb6,0x38,0xef,0xba,0xeb,0xdb,0x31,0xcc,0xd2,0x9b,0xb1,0x83,0xe1,0x15,0x70,0x26,0x6b};
    std::uint8_t product88[32]    = {0x3e,0x4c,0x85,0x55,0xb9,0xbd,0xa4,0x08,0x87,0x66,0x1e,0xcc,0x22,0x1c,0x59,0x6c,0x71,0xdf,0x75,0xd7,0xb6,0x63,0x99,0xa5,0x37,0x63,0x07,0xc2,0x2a,0xe0,0x4c,0xd6};
    mont256_from_bytes(&ctx, multiplier, multiplier88);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product88, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier89[32] = {0xea,0x59,0x67,0x9a,0xed,0x3a,0x32,0xa8,0x6a,0xf2,0x57,0x48,0x8d,0x95,0x9c,0x31,0xfe,0x8a,0xd4,0xa1,0x56,0xd2,0xa6,0x8c,0x02,0xf4,0xb3,0x42,0x74,0x2a,0x80,0x63};
    std::uint8_t product89[32]    = {0xd4,0xb2,0xcf,0x36,0xda,0x74,0x65,0x50,0xd5,0xe4,0xae,0x91,0x1b,0x2b,0x38,0x63,0xfd,0x15,0xa9,0x43,0xad,0xa5,0x4d,0x17,0x05,0xe9,0x66,0x84,0xe8,0x55,0x00,0xc7};
    mont256_from_bytes(&ctx, multiplier, multiplier89);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product89, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier90[32] = {0xf0,0x29,0x05,0x31,0x3d,0x0a,0x27,0x0b,0xb5,0xa4,0x32,0xcf,0x86,0xe3,0xe7,0x26,0x0b,0x0f,0x87,0x3b,0x21,0x14,0xe0,0x68,0x9f,0x27,0xf5,0x2c,0x44,0x92,0x74,0xd2};
    std::uint8_t product90[32]    = {0xe0,0x52,0x0a,0x63,0x7a,0x14,0x4e,0x17,0x6b,0x48,0x65,0x9f,0x0d,0xc7,0xce,0x4c,0x16,0x1f,0x0e,0x77,0x42,0x29,0xc0,0xd0,0x3e,0x4f,0xea,0x58,0x89,0x24,0xe9,0xa5};
    mont256_from_bytes(&ctx, multiplier, multiplier90);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product90, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier91[32] = {0xee,0xa7,0xbb,0x64,0x33,0xa7,0x15,0x68,0x2e,0x5f,0x95,0x0c,0x0c,0xe5,0xaf,0x69,0x43,0x0b,0x91,0xed,0x29,0x54,0xba,0x5c,0xf8,0x1e,0x54,0xdd,0x1c,0x05,0x02,0xc6};
    std::uint8_t product91[32]    = {0xdd,0x4f,0x76,0xc9,0x67,0x4e,0x2a,0xd0,0x5c,0xbf,0x2a,0x18,0x19,0xcb,0x5e,0xd2,0x86,0x17,0x23,0xdb,0x52,0xa9,0x74,0xb8,0xf0,0x3c,0xa9,0xba,0x38,0x0a,0x05,0x8d};
    mont256_from_bytes(&ctx, multiplier, multiplier91);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product91, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier92[32] = {0x72,0x18,0x88,0xff,0x4a,0x3a,0xdf,0x99,0x34,0xb3,0xff,0x60,0xc2,0x6e,0x7a,0x42,0x87,0xf5,0x3d,0xdd,0x4e,0x14,0xd5,0x71,0xa0,0xf0,0x96,0xda,0x4f,0xde,0xbb,0xec};
    std::uint8_t product92[32]    = {0xe4,0x31,0x11,0xfe,0x94,0x75,0xbf,0x32,0x69,0x67,0xfe,0xc1,0x84,0xdc,0xf4,0x85,0x0f,0xea,0x7b,0xba,0x9c,0x29,0xaa,0xe3,0x41,0xe1,0x2d,0xb4,0x9f,0xbd,0x77,0xd8};
    mont256_from_bytes(&ctx, multiplier, multiplier92);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product92, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier93[32] = {0xfe,0x97,0x7c,0x56,0x04,0xa6,0x56,0x51,0xcd,0xbd,0xe7,0x47,0x58,0xd5,0x0f,0x1b,0x45,0x40,0xf4,0x26,0x2d,0x8a,0xd8,0xc0,0xac,0x12,0x7e,0x93,0x80,0x05,0xce,0x74};
    std::uint8_t product93[32]    = {0xfd,0x2e,0xf8,0xad,0x09,0x4c,0xac,0xa3,0x9b,0x7b,0xce,0x8e,0xb1,0xaa,0x1e,0x36,0x8a,0x81,0xe8,0x4d,0x5b,0x15,0xb1,0x80,0x58,0x24,0xfd,0x27,0x00,0x0b,0x9c,0xe9};
    mont256_from_bytes(&ctx, multiplier, multiplier93);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product93, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier94[32] = {0xfa,0x61,0x97,0x74,0x8d,0x11,0x8e,0x37,0x81,0x72,0x8a,0x07,0xbb,0xab,0x27,0xf6,0x04,0xb8,0x15,0x7d,0x03,0xed,0xb9,0x20,0x09,0x75,0x83,0x40,0x40,0x1d,0x68,0xfb};
    std::uint8_t product94[32]    = {0xf4,0xc3,0x2e,0xea,0x1a,0x23,0x1c,0x6f,0x02,0xe5,0x14,0x0f,0x77,0x56,0x4f,0xec,0x09,0x70,0x2a,0xfb,0x07,0xdb,0x72,0x3f,0x12,0xeb,0x06,0x80,0x80,0x3a,0xd1,0xf7};
    mont256_from_bytes(&ctx, multiplier, multiplier94);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product94, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier95[32] = {0xa8,0x87,0xae,0x22,0x1b,0x35,0x41,0x1b,0x72,0x72,0x3b,0x9c,0xef,0x44,0xc0,0xd5,0x3e,0xe4,0xda,0x5a,0x79,0x89,0xe9,0xd0,0x83,0xa4,0xe6,0x29,0x30,0x80,0x38,0x89};
    std::uint8_t product95[32]    = {0x51,0x0f,0x5c,0x45,0x36,0x6a,0x82,0x36,0xe4,0xe4,0x77,0x39,0xde,0x89,0x81,0xaa,0x7d,0xc9,0xb4,0xb5,0xf3,0x13,0xd3,0xa0,0x07,0x49,0xcc,0x52,0x61,0x00,0x71,0x13};
    mont256_from_bytes(&ctx, multiplier, multiplier95);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product95, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier96[32] = {0xe3,0x83,0x8b,0x9e,0xd5,0xa9,0x42,0x2a,0x8b,0xc0,0x83,0x11,0x7e,0xb8,0x6c,0x57,0xa8,0x11,0x00,0xa1,0x6e,0xa3,0x30,0xa1,0xa6,0x6d,0x58,0xb5,0xd1,0xa4,0xc0,0x1e};
    std::uint8_t product96[32]    = {0xc7,0x07,0x17,0x3e,0xab,0x52,0x84,0x55,0x17,0x81,0x06,0x22,0xfd,0x70,0xd8,0xaf,0x50,0x22,0x01,0x43,0xdd,0x46,0x61,0x42,0x4c,0xda,0xb1,0x6b,0xa3,0x49,0x80,0x3d};
    mont256_from_bytes(&ctx, multiplier, multiplier96);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product96, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier97[32] = {0x3a,0xc4,0xda,0x9a,0xfb,0x81,0x39,0x21,0x37,0x16,0x1c,0x16,0xb0,0x0f,0xd7,0xbb,0x4e,0xca,0xde,0xa2,0x81,0xb6,0x2b,0xb5,0xf8,0x66,0x64,0xae,0x64,0xa1,0x49,0xf5};
    std::uint8_t product97[32]    = {0x75,0x89,0xb5,0x35,0xf7,0x02,0x72,0x42,0x6e,0x2c,0x38,0x2d,0x60,0x1f,0xaf,0x76,0x9d,0x95,0xbd,0x45,0x03,0x6c,0x57,0x6b,0xf0,0xcc,0xc9,0x5c,0xc9,0x42,0x93,0xea};
    mont256_from_bytes(&ctx, multiplier, multiplier97);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product97, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier98[32] = {0x23,0xc4,0x9c,0xae,0xa2,0xcf,0x62,0xba,0xba,0x95,0x88,0x10,0xb4,0xeb,0xf4,0xb6,0xe1,0xc6,0x0a,0xa3,0xd5,0x10,0xbb,0x04,0x32,0xd9,0x0d,0xcd,0x57,0xbb,0x7d,0x97};
    std::uint8_t product98[32]    = {0x47,0x89,0x39,0x5d,0x45,0x9e,0xc5,0x75,0x75,0x2b,0x10,0x21,0x69,0xd7,0xe9,0x6d,0xc3,0x8c,0x15,0x47,0xaa,0x21,0x76,0x08,0x65,0xb2,0x1b,0x9a,0xaf,0x76,0xfb,0x2e};
    mont256_from_bytes(&ctx, multiplier, multiplier98);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product98, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
    std::uint8_t multiplier99[32] = {0x03,0xa6,0x39,0x66,0x21,0x3b,0xca,0x7f,0xd6,0x44,0xde,0x2f,0x0d,0xec,0x68,0x23,0xfb,0x5c,0x9d,0x56,0x58,0xf9,0x2d,0xea,0xfd,0x4b,0xd0,0x30,0x67,0x9a,0x44,0xdd};
    std::uint8_t product99[32]    = {0x07,0x4c,0x72,0xcc,0x42,0x77,0x94,0xff,0xac,0x89,0xbc,0x5e,0x1b,0xd8,0xd0,0x47,0xf6,0xb9,0x3a,0xac,0xb1,0xf2,0x5b,0xd5,0xfa,0x97,0xa0,0x60,0xcf,0x34,0x89,0xba};
    mont256_from_bytes(&ctx, multiplier, multiplier99);
    mont256_dbl(&ctx, product, multiplier);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product99, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_dbl");
    }
}
